Dynamical conductivity at the dirty superconductor-metal quantum phase transition 
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We study the transport properties of ultrathin disordered nanowires in the neighborhood of the 
superconductor-metal quantum phase transition. To this end we combine numerical calculations 
with analytical strong-disorder renormalization group results. The quantum critical conductivity 
at zero temperature diverges logarithmically as a function of frequency. In the metallic phase, it 
obeys activated scaling associated with an infinite-randomness quantum critical point. We extend 
the scaling theory to higher dimensions and discuss implications for experiments. 



Electrical transport in low-dimensional strongly fluc- 
tuating superconductors has been the subject of intense 
experimental investigation for almost half a century [1]. 
Recently, advances in experimental techniques have al- 
lowed for the fabrication of ultrathin metallic nanowires 
having diameters smaller than the bulk superconduct- 
ing coherence length, but large enough to include many 
transverse channels for electronic conduction. Resistance 
measurements have shown that the thicker among these 
wires exhibit a well-defined phase transition from a re- 
sistive to a superconducting state with decreasing tem- 
perature, while thinner wires appear to remain resistive 
down to the lowest temperatures measured [2-6]. 

It has been proposed [7-10] that these experiments 
may be understood in terms of a superconductor-metal 
quantum phase transition (SMT) driven by pair-breaking 
interactions, possibly due to random magnetic moments 
trapped on the wire surface ]11]. A description of this 
transition is provided by a theory, first proposed by 
Feigel'man and Larkin ]12], of a complex Cooper pair 
order parameter whose fluctuations are damped by de- 
cay into unpaired electrons ]13-17]. 

As the nanowires are prone to random variations in di- 
ameter and because of the random positions of the pair- 
breaking moments, quenched disorder plays an impor- 
tant role. The thermodynamics of the disordered SMT 
has been analyzed both analytically ]18] and numerically 
]10] in the relevant case of one space dimension. It is 
governed, for any nonzero disorder strength, by a non- 
perturbative infinite-randomness critical point (IRCP). 
This IRCP is in the same universality class as the mag- 
netic quantum critical point of the random transverse- 
field Ising chain despite the fact that the two systems 
have different symmetries: The clean transverse-field 
Ising chain can be described by relativistic free fermions 
(and, therefore, dynamical exponent z = 1) whereas the 
clean SMT is described by overdamped 0(2) fluctuations 
with z = 2. The homology lies in the marginal dynamics 
of finite size clusters in both models [20] which are the 
famous rare regions of Griffiths-McCoy physics [21]. 



Many asymptotically exact results for the random 
transverse-field Ising chain [19, 22[ apply directly to the 
SMT via universality. The IRCP is characterized by ac- 
tivated dynamical scaling: Ln ~ [ln(rio/r2)]^/'''. Here, H. 
is the characteristic energy of the order-parameter fluc- 
tuations on length scale Ln, is a high-energy reference 
scale, and = 1/2 is known as the tunneling exponent. 
The exponential length-energy relation implies that the 
dynamical exponent z is formally inflnite. Moreover, 
the magnitude of the order-parameter fluctuations fi also 
scales logarithmically with energy, ^ [ln(J7o/f^)]'^, 
where the cluster exponent (f) = {1 + \/5)/2 is the golden 
ratio. Approaching criticality, the correlation length di- 
verges as ^ ~ where i/ = 2 and 6 measures the 
relative distance to the critical point. 

In this Letter, we study experimentally important 
transport properties at the pair-breaking SMT of disor- 
dered nanowires. We report both analytical and numer- 
ical calculations of the zero-temperature finite-frequency 
Aslamazov-Larkin [23[ fluctuation corrections to the con- 
ductivity a (w). At criticality, the real part of the conduc- 
tivity diverges as (T'(a;) ^ [\n{uJo/u})]^^^' with vanishing 
frequency co {loq is a reference frequency). Off criticality, 
it satisfles the unconventional activated scaling form 



where ^^ix) is a universal scaling function. In the re- 
mainder of this Letter, we sketch the derivation of these 
results and discuss their experimental implications. 

We begin by introducing a one-dimensional continuum 
model of Cooper pairs in the presence of Ohmic dissipa- 
tion and disorder at T = [7, 8, 10, 18[ 



S = dx dr 



D{x)\d,.'<Sj{x,T)\^ + aix)\^ix,T)\'' 



+ ^|*(x,r)|^] +/dxy ^j{x)\oj\\^{x,u)\\ (2) 

Here, '^{x,uj) is the Fourier transform of 4'(x,t), a com- 
plex superconducting order parameter at position x and 
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imaginary time r. We have explicitly included the ran- 
dom spatial dependence of all coupling constants. Sta- 
bility and causality constrain u,^{x) > 0, and we may 
choose a gauge where D{x) > 0. The quantum phase 
transition is tuned via S = a — ac [10]. 

To proceed, we use a lattice discretization of the con- 
tinuum action (2) in the limit of a large number of order- 
parameter components. This limit has no impact on 
the character of the critical point [10, 18] and leads to 
a quadratic action 

'5° = E/ + mj)^,iu:), (3) 

hi 

where the coupling matrix My = {Dj/ ^ji'yj)Afj + 
{T'iHi)5ij (with the discrete nearest neighbor Lapla- 
cian) must be determined self-consistently by solving 
n=a, + and ^i{uj) ^ has 

been rescaled. Note that the full effects of disorder can 
be realized while fixing 7^ = 7 to be constant [24[. 

We are now in a position to directly evaluate the dy- 
namical conductivity via the Kubo formula [25[ 



10' 
10" 

10-1 

10-5 

10-7 



128 



t 



1 



• 0.93 - 
■ 0.68 : 

• 0.43 - 

▲ 0.28 - 

▼ 0.18 - 



, ^ t, ^ f Q 

10-*^ 10-^^ 10--* 10-^ 10-- 10-' 10" 10' 10^ 10-^ lO** 



Figure 1: (Color online) The disorder averaged real conduc- 
tivity for chains of 128 sites as a function of frequency mea- 
sured in terms of a UV cutoff ojq for different values of the 
distance from the critical point, 5. 
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where the current is given by Jj{T) = 
(2zeM)Z?, [vl/*(r)vI/^.+,(r)-vI/*^i(r)vI/^.(r)] with dia- 
magnetic contribution V = (8eV7^) A (l^'i(O)P). 
The Kubo formula can be evaluated by employing the 
spectral decomposition of M in terms of its eigenvector 
Vij and accompanying diagonal eigenvalue = eiSij 
matrices defined by Et ^ik^kj = 
real part of the conductivity reads 



The resulting 



where 

ICa.bi^) = 



~h 



a, 6 i,j 



(4) 



{-2ea(e^ 



- o;^) arctan(aj/eQ) 



- 2eb{-el + 4+ arctan(cj/eb) 



(5) 



The validity of the scaling form (1) can be tested via 
a full numerical evaluation of Eqs. (4) and (5). This is 
possible by exploiting an efficient algorithm for comput- 
ing the self-consistent pairing eigenmodes of So for large 
system sizes [10] . We have evaluated the conductivity (4) 
for chains of various length with up to 128 sites averaged 
over 3000 disorder realizations. For clarity, we limit our 



analysis to the largest size, L ~ 128, as the extrapola- 
tion to the thermodynamic limit L — )■ (X) is non-trivial 
due to the crossover between ^ and L for the range of i5 
considered here. The results are displayed in Fig. 1. 

For probe frequencies uj much larger than the char- 
acteristic fluctuation energy scale wq of the chain, we 
fully saturate all quantum dynamics and observe a triv- 
ial a'{uj) ^ 4e^/(/iw^) conductivity. On the other hand 
for a; ^ Wo, cr' appears to be suppressed by a (S-dependent 
exponent. As we approach criticality {S — > 0) the func- 
tional form of the average conductivity is not easily as- 
certained from Eq. (4) due to the softening of critical 
modes, but the apparent disappearance of this exponent 
is fully consistent with the scaling theory in Eq. (1). 

The predictions of Eq. (1) can be more thoroughly 
confirmed by searching for consistent scaling of the data 
shown in Fig. 1 after dividing by [ln(a;o/w)]^/'^ and re- 
plotting as a function of the dimensionless scaling vari- 
able X = S^'^ \n{LjQ/Ljj). We find excellent data collapse 
over five orders of magnitude as shown in Fig. 2. 

From the scaled data, we extract the universal pref- 
actor of the conductivity at the critical point, $(j(0) ~ 
0.70(4). Furthermore, an empirical analysis of the nu- 
merical scaling function in Fig. 2 suggests that ^cr{x — 
00) ^ x~'^c~-^^ where A ^ 0(1). The relevant limits 
of $(a;) can also be inferred by appealing to the naive 
scaling prediction that in d = 1 the conductivity should 
be equal to Ae^/h multiplied by a length. Activated dy- 
namical scaling dictates that at criticality, lengths scale 
like [ln(wo/'^)]"'^^''^ whereas in the Grifhths phase the rel- 
evant length scale is the correlation length ^ ~ I'^l^'', 
fixing the x-dependent power in front of the exponential 
in <f>cr(a; — > 00) in order to cancel the logarithmic prefac- 
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Figure 2: (Color online) The disorder averaged real conduc- 
tivity scaling function $ct as a function of the dimensionless 
scaling variable x = S'^'^ ln{ujo/uj) for different values of 5 as 
criticality is approached from above. The line is a guide to 
the eye showing the probable functional form of $o- (see text). 



tor in Eq. (1). 

Let us now turn to an analytical derivation of the dy- 
namical conductivity. Near an IRCP, the conductivity 
will be dominated by large rare regions which are locally 
in the superconducting phase, i.e., by small clusters with 
exceptionally strong links D and typically small gaps r. 
In the low frequency limit, the effective links and gaps 
of these clusters can be quantified by a renormalization 
group analysis [18] that we invoke later. For now, we 
approximate each dominant cluster as a single two-site 
system with ri.2 being the effective local gaps and D the 
effective link strength. In such a simple model, the con- 
ductivity Eq. (4) can be evaluated exactly, resulting in 



8eW 

'■2sitc — 2 '^12lWj, 



(6) 



where the eigenvalues of the 2x2 coupling matrix are 

2ei,2 = [D/-f + in + r2)/7 ± VWW+Tn^^WFF] 
and we have introduced £ as the length of the link con- 
necting the two sites. Due to the presence of the fac- 
tor D"^ in Eq. (6), the average over the contributions of 
all two-site clusters will be dominated by those with an 
anomalously large links D ^ ri.2. Hence, we may evalu- 
ate the conductivity by averaging 



a' (fl) 



dD 



1 pD 

deP{D,£) dri?(r)aLito, 

^0 



(7) 

where ft = joj is the energy scale at which effective clus- 
ters with gap r connected by links of magnitude D and 
length i appear with probability density R{r)P{D,£), 
and HQ is the density of such clusters. The final step 
consists of using the asymptotic value of /Ci2(w) ^ 



uj ^ln[r2/r]. in the strong disorder limit and i^,D ^ r 
where r = (ri + r2)/2. Using the values of P{D,£), R{r) 
and nn at criticality [19, 22[: P{y) = Jd£P{y,£) = 
e-y, with y = \n{n/D)/\n{nQ/n), R{y) = P{y) as 
well as the relation n^^ ~ {£) - / ' be- 

tween density and average separation of clusters, we ar- 
rive at a' - nn{{D/nf£^) (ln(17/r)) - [\n{nn/n)Y''l' , 
where = 2 recovering Eq. (1) at criticality where 
5""^ \n{VLQ/U) — > 0. Correlations between the links D 
and their lengths £ have no effect on the leading logarith- 
mic divergence of this result at criticality. The analysis 
of Eq. (7) in the metallic and superconducting Griffiths 
phases crucially depends on the careful treatment of the 
correlations between £ and D. The resulting expressions 
are quite involved and will be discussed elsewhere [26[. 

We now discuss a sub-leading correction to the scal- 
ing of cr' {lo). In Eq. (7), the dissipative z — 2 dy- 
namics causes the relation between energy and the 
measured frequency u to have a logarithmic correction 
[18[ = 7o/io<^, where 70 is the bare dissipative cou- 
pling and /io is the mean value of superconducting order- 
parameter fluctuations. At criticality /in ~ [ln(i7o/i^)]'^ 
and ln(f2o/i^) = ln(wo/i^) up to log(log) corrections. 
In the metallic Griffiths phase, where presumably any 
real experiments on metallic nanowires would take place, 
/^a ^ J"'^'-^"'^-' ln(f2o/^^) and the logarithms of energy 
and frequency are no longer simply equivalent. However, 
the exact value of ^n can be obtained from the imagi- 
nary part of the dynamical order-parameter susceptibil- 
ity [10, 18[, and its inclusion leads to quantitatively better 
data collapse as we extend the scaling theory of Eq. (1) 
deeper into the Griffiths regime [26[. 

In order to place the appearance of the link length £ 
in Eq. (6) and the average of Eq. (7) on firmer footing, 
we invoke the real space renormalization group technique 
of Refs. [18, 27], providing direct access to the renormal- 
ization of the current operator. Starting with a chain 
described by the effective action of Eq. (3), we proceed 
by searching for the largest local coupling of the chain, 
17 = max{ri,I?i}. Suppose (i) VL = r2, site 2 is then 
strongly fluctuating and can be integrated out of the the 
system leading to an effective coupling D = DiD^jr^ 
between sites 1 and 3. Their relative distance is given by 
£ = ^1 +^2, and the total current through clusters 1 — 3 is 
£iJi+£2J2 = where J = [2ie/-ih)b[^l'^r^-'^l^{\. If 
on the other hand (ii) Q. = D2, sites 2 and 3 are strongly 
coupled forming an effective cluster 2 in which the ef- 
fective gap is f2 = T2f3/D2. The total current is then 
iiJi + ^2^2 + 4^3 = hJi + £2J2, where £1 = £1 + i^, 
£2 = 1^2 + ^3, k = {2iehh)Di[^l-^^^ - and 
J2 = (2ie/7/i)D3[*?*4 - vl'l*^]- After process (i) or 
(ii), the energy scale O is lowered and the disorder of 
the effective system is increased [19]. Iterating this pro- 
cedure leads to the probability distribution of gaps and 
links connecting the effective clusters as a function of en- 
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ergy [18] providing formal justification of Eq. (7). 

We now compare our results with the transport prop- 
erties of other systems governed by infinite randomness. 
To the best of our knowledge, only one other such sys- 
tem has been studied: spin conductivity in the dimer- 
ized antiferromagnetic spin-1/2 chain [27]. At criticality, 
the spin conductivity also diverges logarithmically, but 
with a weaker power than here, and is found to obey the 
scaling form crspin(w) ~ ln(a;o/w)$spin ['51n(cJo/<^)], where 
$spin(a; -)> 0) ft! 7/180 and $spin(a; oo) cce"^"". Al- 
though the thermodynamics of our system and the spin 
chain are funneled into the same universality class by 
the disorder, their transport properties are not universal 
because their underlying dynamics are different. 

Finally, we highlight that the methods discussed here 
should also apply in higher dimensions. Specifically, at 
criticality, the system will also be governed by an infinite- 
randomness critical point [18, 22] but with different expo- 
nents "0, V, and (/). The dynamical conductivity will like- 
wise be dominated by rare and locally superconducting 
strongly coupled regions. Evaluating Eq. (7) immediately 
leads to (t'(w) ^ [\n{ujo/uj)]''^~''-^^'^' , since the spatial di- 
mension enters explicitly only via the density of clusters 
~ L^"^, and the average conductivity of a two-site 
cluster (cTjsit;,) ^ [ln(wo/'^)]^^'''- As an immediate con- 
sequence, in the limit w — > 0, the critical conductivity 
vanishes for d > 2, and becomes constant at d = 2. 

In conclusion, we have presented the numerically eval- 
uated dynamical conductivity cr'(a;) for a model believed 
to describe the physics of disordered nanowires close to 
a superconductor-metal quantum phase transition and 
placed it in an analytical framework computed via the 
strong disorder renormalization group. We have shown 
that (j'{uj) diverges logarithmically as |lna;|^ at critical- 
ity and obeys scaling in the metallic (Griffiths) phase 
with asymptotics dictated by naive dimensional analy- 
sis of physical quantities. Our results may be directly 
applicable to experimental transport measurements on 
thin nanowires which remain metallic as T — > 0, such as 
those reported in Ref. [6]. By studying wires of vary- 
ing thickness at low temperatures, it may be possible 
to reach the critical regime where the logarithmic diver- 
gence of the fluctuation correction to the conductivity for 
uj > T could be directly observed. In d = 2, we expect 
a' to be frequency- independent, consistent with experi- 
ments in disordered thin superconducting films (see, e.g. 
Ref. [28]). Finally, by presenting solid predictions arising 
from an effective action of strongly repulsive dissipative 
Cooperons, we lay open an avenue for the experimental 
investigation of the efficacy of such models when applied 
to dirty low dimensional superconductors. 
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